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1. Introduction 
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Viscous GCG(generalized Chaplygin gas) cosmology is discussed, assuming that there is bulk vis- 
cosity in the linear barotropic fluid and GCG. The dynamical analysis indicates that the phase 
Wg — — 1 + V3'yKTg/{'y — \/3ktj) is a dynamical attractor and the equation of state of GCG 
approaches it from either Wg > — 1 or Wg < —1 depending on the choice of its initial cos- 
mic density parameter and the ratio of pressure to critical energy density. Obviously, the equa- 
tion of state Wg can cross the boundary Wg = —1. Also, from the point of view of dynam- 
ICS, the parameters of viscous GCG should be in the range of 7 > \/3K,r-y/(l — \/3KTg) and 
O. 0<.<l + 73.W(7-^/3...-V3....). 

(N 

CMB anisotropy, supernovae and galaxies clustering strongly indicate that our universe is spatially flat, with two 
CNJ ■ thirds of the energy content consisting of dark energy, a substance with negative pressure which can make the universe 
I expand in an acceleratirig fashion. Present observation data constrain the range of the equation of state of dark energy 
. as —1.38 < w < —0.82 which indicates the possibility of dark energy with w < — 1. If this is so, the universe can 
have some strange properties such as a future finite singularity which has been dubbed Big Rip 2] or Big^niashQ. 
Proposed candidates for dark energy include the cosmological constant, quintessence with a single fieldQ or with 
N coupled fields[3, phantom field with canonical^ or Born-Infeld type Lagrangian0, A;-essence|^ and generalized 
IT) I Chaplygin gas(GCG)[3 which is based on the Chaplygin gas0. In particular, we have extended the equation of 
state of GCG to the regime w < —1 regime and shown that the GCG parameter is constrained by the dynamics to 
lie in the range < a < 1 

O I On the other hand, the role of dissipative processes has been conscientiously studied^3|- Dissipative effects, 
I including both bulk and shear viscosity, play a very important role in the evolution of the universe. The viscosity 
theory of relativistic fluids was flrst suggested by Eckart^lSj and Landau and Lifshitz|ljj, who considered only flrst- 
^ order deviation from equilibrium, which lead to parabolic differential equations and hence to an infinite speed of 
, propagation for heat fiow and viscosity, in contradiction with the principle of causality. The relativistic second-order 
IL^ theory was founded by Israelii^ and developed by Israel and Stewart [l6|. and has also been used in the evolution of 
] the early universe|l7j|. However, the character of the evolution equation is very complicated in the framework of the 
. full causal theory. Therefore, the conventional theory|Q| is still applied to phenomena which are quasi-stationary, 
?H i.e., slowly varying on space and time scales characterized by the mean free path and the mean collision time. In the 
case of isotropic and homogeneous cosmologies, the dissipative process can be modelled as a bulk viscosity C within a 
thermodynamical approach. The shear viscosity 77 will be neglected, which is consistent with the usual practice 18j. 
The bulk viscosity introduces dissipation by only redefining the effective pressure, Peff, according to Peff — P — 3^i/ 
where C is the bulk viscosity coefficient and H is the Hubble parameter. The condition C > assures a positive entropy 
production in conformity with the second law of thermodynamics |0. We are interested in the case 
where r is a constant. This assumption implicates that C is directly proportional to the divergence of the cosmic 
fluid's velocity vector. Therefore, it is physically natural, and has been considered previously in an astrophysical 
context (23. 

In the present paper, we consider a viscous GCG cosmological model for the expanding universe, assuming that 
there is bulk viscosity in the linear barotropic fluid and GCG. The dynamical analysis indicates that the phase 
Wg = —1 + ^/3jKTg/{j — ^/SkTj) is a dynamical attractor and the equation of state of GCG approaches it from either 
Wg > —1 or Wg < ~1 depending on the choice of its initial cosmic density parameter and the ratio of pressure to 
critical energy density, where we assume y/SnTg and y/SnTj are small compared to 1 . If indeed dark energy with 
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ui < — 1 is within the regime of possibilities, then it would seem inevitable to inquire about a transition from w < —1 
to w > —1, although a generalized scalar field cannot safely traverse w = —1. It is clear that w can cross the boundary 
Wg = ~1 in the viscous GCG cosmology. Also, we have used cosmological dynamics to place constraints on some of 
the parameters. 

2. Viscous GCG 

GCG has a very simple equation of state, pg — — Af^^^+^^/p^ , which yields an analytically solvable cosmological 
dynamics if the universe is GCG dominated. Another advanta ge o f introducing GCG is to unify dark energy and dark 
matter into one equation of state, also known as a uartessence |2 ij . However, detailed numerical analysis turns out to 
disfavor the dark matter modelled by the GCG equation of state |23. But no observation has so far ruled out the 
possibility of GCG as dark energy. Therefore, it is quite possible that our universe contains a dark energy component 
modelled by the GCG as well as another linear barotropic fluid component with the equation of state p-y = (7 — 
However, in this section, we focus first on the viscous GCG system, as the cosmological dynamics are analytically 
solvable if GCG is dominant . In the flat FRW universe, the field equations are 

= = ^Pg (1) 

2 

H + H^ ^ ---^(Ps + 3pg-9iJC) (2) 



a 6 



The conservation equation is 



Pg+3H{pg+Pg-3HC)=0 (3) 

Using the GCG equation of state, we obtain 



a dpg M4("^ 



3 da ' 

We shall be interested in the evolution of the late universe, from t = to onwards. From Eq.(4), we have 



-V3kp]C = (4) 



a = ao \ exp 



(5) 

II I i-H I t -1- I J I ./■_/>,. ^l-i I TT /( ^ 1 ^( y -I- I / 



which is the general relation between the cosmological scale factor a and the energy density pg. 

In what follows we investigate two different solvable cases. If we choose <^ = Tgy/p, the energy density is given by 



Pa 



KTn 



Obviously, | pg |— > cx3 or will not happen when a ag ^ 0, corresponding to no future singularities. As a — > oo, we 
have Wg —(1 — ^/inTg) > —1. 

As a second case, if we take C = Tgp"'^'^ then Eq.(4) reads 

_^^(,;«,_,j«+M-«. + V3»...-«. (7) 

and if Tg < (4\/3KAf'*("+^))^^, it has two stationary points at pg — ( ^.-^ ) ^ = pi and pg — ( ) ^ = P2 

where A = (1 — 4\/3KTgAi^''("+^')5. The first stationary point is an attractor, and corresponds to a ^ oo. The second 
stationary point is unstable, and corresponds to a = 0. If po < P\ or pi < po < P2 then the model evolves As a — > oo, 

Pg ^ Pi and Wg r— A ^ However, it po > P2 then as a ^ aol "+i „+i I , p„ ^ oo and 
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— > cxD so the model encounters a Type III singularity (where the pressure and density diverge in finite proper time 
but the metric remains non-singular). 

3. Autonomous system 

A general study of the phase space system of quintessence and phantom in FRW universe has been given in 
R.efs. |2^l23 |. For the viscous GCG cosmological dynamical system, the corresponding equation of motion and Einstein 
equations can be written as 



H = - — {p^+p^^SC^H + pg+pg-SCgH) (8) 
= -3H{p.^+p^-3C-tH) (9) 

Pg = -3HiPg+Pg-3CgH) (lO) 

= y(P7 + Pff) (11) 

where pj is the density of fluid with a barotropic equation of state p-y = {j — ^)p-y, and < 7 < 2 is a constant that 
relates to the equation of state by t/;^ = 7 — 1; pg and Pg are the energy density and pressure of GCG, respectively. 
The overdot represents the derivative with respect to cosmic time t. We choose Cg — '''gy/P — VSn^^TgH and 
— T^Y^ = \/in~^T^H . 

To analyze the dynamical system, we rewrite the equations with the following dimensionless variables: 

" 3i?2 

N = Ina (12) 



The dynamical system can be reduced to 



dN 

From Eq.(ll), we have 



= —i{x + y) + 2>V?tKTg + ix['y{l — x) + X + y\ — iVii^iTg + Ty)x (13) 
= ia{y H ) — iVSaKTg — h 3y[7(l — x) + x + y] — 3\/3K(Tg + T^)y (14) 



ng + n^ = i (15) 



where fig = x and fl.y = 3^3" the cosmic density parameters for GCG and linear barotropic fluid, respectively. 
The equation of state can be expressed in terms of the new variables as 



and the sound speed is 



Wg = — ^- (16) 
Pg X 



cl = -a^ (17) 

X 



In the following we need to consider the two cases of 7 = 1 and 7 7^ 1. In the case of 7 7^ 1, the critical points of 
the system are 
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2(7-1) 

and 

V^KT^ _ /- \/3i 

^,-l + V3KTg + 

7 7 



(.(^\,(^)) = (l - ^,-1 + V3«r, + ^) (20) 



which correspond to the Uncar barotropic fluid-dominated phase, GCG matter-dominated phase and GCG vacuum- 
energy-dominated phase, respectively. Here 

S = [4V3(7 - l)KTg + (7 - 1 - VSKTg - VSnT^f]^ (21) 

If we linearize the system near the critical points (a;'-*\ y^*^), i = 1, 2, 3 and then translate the system to the origin, we 
can readily write the first order perturbation equation as 

where C/ is a 2-column vector consisting of the perturbations of x and y. A'*' is a 2 x 2 matrix for the critical point 
(x^*\ 2/^*^). The stability of the critical points is determined by the eigenvalues of the matrix A^'^ at the critical point 
(a;(*\y(*)). For the point {x^^\y^^^), the two eigenvalues are 

A^^^ = 3S (23) 
^(1) _ 3(1 + a){j - 1 + V^KTg - 2\/^^KTg - a/Skt^ - S) 
^ 7 — 1 — \/3KTg — \/3ktj — S 

Obviously, the linear barotropic fluid-dominated phase is unstable, and so evolves to GCG-dominated phase. For the 
point {x^'^\y^^^) the two eigenvalues are 

Af^ = -3S (24) 

,(2) _ S{1 + a){'-f — 1 + V^KTg ~ 2^^KTg — VSkt^ + T,) 
Ao 



^2 



7- 



1 — \/3KTg — VSKT-y + S 



It can be shown that A22) > 0, and so the point is stable if a > — 1 and 7 < ^ ■ However, since we want the 

GCG to first behave as matter and then evolve to behave like dark energy, it will not be appropriate if (x^^^y^"^^) 
corresponds to a stable attractor phase. In other words 7 < fl'^"'' should not be considered in the real models. For 
the critical point (a;^^\ t/^^^), the corresponding eigenvalues of matrix AS"^^ are 

\f = %= (e + S) (25) 

2(7 — V^KT-y) 

Af = — (e-s) 

2(7 - VSkt^) 



where 



e = 7^ + [(1 + a)(l - VSKTg) - 2\/3kt^]7 - (1 + a)V^KT^ + SK^iTgT-, + r^) (26) 

1 

+ [7^ + [(1 + a){l- V^KTg) - 2V3ktj]^ - (1 + a)\/3KT-y + 3K^{TgT-y + T^)] ^) ' 
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TABLE I: The initial values of x and y in the plots FIG.1-FIG.3 
X 0.14 0.15 0.16 0.17 0.18 0.18 



y -0.19 -0.18 -0.17 -0.16 -0.15 -0.14 

TABLE II: The initial values of x and y in the plot FIG.4 

ITOXiim 0.0240 0.0340 0.0440 0.0540~ 
y -0.0042 -0.0036 -0.0030 -0.0024 -0.0018 

Note that A^^ Aj^^ = 9(a -|- 1)(7 — V^jKTg — V^kt-j), and it can be shown that (a;^^\ y^^') is stable for a > —I and 

y > . 

In the case of 7 = 1, the system has only two critical points, the first one is {x^*\y^*^) = (^ rg+T^ ' ^) ' second one 

is still Eq.(20) with everywhere replaced by 7 = 1 . At the first critical point, the eigenvalues are A^*^ = —3V3K{Tg+T-y) 

and Aj*' = 'S{a + 1)(1 — \/3K{Tg + r^), and so the point is unstable if a > — 1 and 1 — V3K{Tg + T-y) > 0. It can be 
shown that the second critical point is stable under the same assumptions. So we obtain the same stable point in 
both cases. This critical point corresponds to a GCG-dominated phase 

_V3kx, (27) 
7 7 

and its equation of state is 

-« = -l + ^%^ (28) 

In the Tj = Tg = case, we have fig = 1 and Wg = —1 which is a late time de Sitter attractor[21]. In the r-y = and 
Tg case, this analysis is consistent with the results of exact solution (6). 

4. Numerical analysis 

Next, wc study the above dynamical system numerically. For dcfinitencss, we choose the parameters to be 7 = 1 and 
a = 0.5. The initial values of x and y are chosen as shown in TABLE I and the results are contained in FIG.1-FIG.3. 
Prom FIG.l, we can observe that all orbits tend to an attractor which corresponds to Wg = — From FIG. 2, we 
can observe that for different initial values of pg and pg, the equation of state Wg approaches the attractor Wg = — |i 
from either Wg > —1 ov Wg < —1. Irrespective of whether the initial choice has Wg > —1 01 Wg < —1, the equation 
of state will eventually mimic that of quintessence in the viscous cosmology. It is worth noting that if we choose the 
parameters so that the GCG behaves as phantom, it is no longer possible to make it behave as matter at an early 
epoch unless the viscosity coefficients are anomalously large. However, in our setup of this paper, we have included a 
linear barotropic fluid that could be used to mimic the matter sector of our universe and thus GCG can be considered 
only as dark energy. The evolution of the sound speed Cg is shown in FIG. 3, where tends to ^ for different initial 
values of pg and Pg, and VSkt^ = 0.03, V^KTg = 0.06. 

The initial values of x and y are chosen as shown in TABLE II and the results are contained in FIG.4. The different 
values of the parameters and Tg are chosen as shown in TABLE III and the results are contained in FIG. 5. 



TABLE III: Tlio parameter values of r,, and r- 
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FIG. 1: The phase diagram of the viscous GCG system in terms of x and y for different initial values of x and y, and 
VSnTg = 0.06, VSkt^ = 0.03. 




FIG. 2: The evolution of the equation of state of viscous GCG for different initial values of x and y, where we have taken 
x/SKTg = 0.06, y/3KT^, = 0.03. The curves from bottom to top correspond to the initial conditions specified in Table I from left 
to right respectively. 



5. Conclusions 



In the present paper, we have analyzed the dynamical evolution of viscous GCG for different parameters and initial 
conditions. Our results are as follows: 

(1) We have shown that different initial values of pg and Pg will lead to different tracks {wg > — 1 and Wg < —1 

) for the equation of state Wg as it approaches the dynamical attractor phase = — 1 H — . That is to say, 

the comparison of theoretical values to observational ones places a constraint on the bulk viscosity coefficient. Recent 
astrophysical data indicate that the effective equation of state parameter Weff lies in the interval: —1.38 < w^ff < 
—0.82 0, so that we have the constraint V^K{^Tg + r^) < 1. 

(2) We can also use the requirement that < < 1 in the GCG to place the constraint < a < IH v ^^'^3 

on the parameter a in the model . 

(3) The equation of state Wg can cross the boundary Wg = —1. 

(4) From the point of view of the dynamics, the bulk viscosity coefficient should satisfy 7 > . 

(5) In the viscous model, the ratio of cosmic density parameters ^ approaches the constant If we suppose 



7 



0.65 




1 2 3 4 5 

N 

FIG. 3: The evolution of the sound speed for different initial values of x and j/ in a viscous cosmological model with 
\/3firg = 0.06. v^KT-y = 0.03. The curves from top to bottom correspond to the initial conditions specified in Table I from left 
to right respectively. 
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FIG. 4: The evolution of the cosmic density parameter for matter flj and fig respectively at different initial values of x and y, 
andVSKTg = 0.02, V^KT-y = 0.01. 



the present universe is in the epoch of approximating to a constant, we can obtain the constraints \/3kt^ = 0.87 
and ^/SnTg < 0.7. 

Finally, it is worth noting that as a phenomenological model of the evolution of the late universe, it is reasonable 
for our model to leave the causal viscosity theory out of account. 
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FIG. 5: The evolution of the equation of state of viscous GCG for different bulk viscosity coefficients, where we have taken 
a; = 0.15 and y = —0.18. The curves from bottom to top correspond to the parameter values specified in Table 3 from left to 
right, respectively. 
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